Solar Physics 
DOI: 10.1007/' 



Oscillations in the Solar Faculae. III. The Phase 
Relations between Chromospheric and Photospheric 
LOS Velocities 

N.I.Kobanovi • A.S. Kustov^ • S.A. Chupini • 
V.A. Pulyaev^ 

© Springer •••• 

o 

Abstract An analysis of line-of-sight velocity oscillation in nine solar faculae was 
O ' undertaken with the aim of studying of phase relations between chromosphere 

I (He I 10 830 A line) and photosphere (Sii 10 827 A line) five-minute oscillations. 

We found that time lag of the chromospheric signal relative to photospheric one 
varies from -12 to 100 seconds and is about 50 seconds on average. We assume that 
the small observed lag can have three possible explanations: i) convergence of for- 
mation levels of Hei 10 830 A and Sii 10 827 A in faculae; ii) significant increase 
■ of five-minute oscillation propagation velocity above faculae; in) simultaneous 

(~| I presence of standing and travelling waves. 
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I 1. Introduction 

Solar faculae, usually called "plages" in the photosphere or calcium and hydrogen 
^ ■ flocculi in the chromosphere, are often observed in the solar atmosphere. They 

00 I may play a significant role in energy exchange processes between layers of the 

. solar atmosphere due to their prevalence. One of possible ways to transport 

energy is waves propagating in faculae. The facular oscillations have been actively 
stud ied since the 1960s (|Orrall, 1965 Howard, 1967_She eley and Bhatnagar, 



1971: 



Deubner, 1974[ |Teske, 1974[ |Cram and Woods, 1980D 



lOrrall, (,1965j) found no evidence of a difference between the velocity ampli- 
tudes underlying weak plages or the chromospheric network, and regions free 
of calcium flocculi. There was only a slight indication that the periods are a 
. ^ ^ little longer in photospheric regions that underlie the network. iHowardi p967p 

■ found that the amplitudes of five-minute oscillations are about 25% weaker in 

H I regions with a significant magnetic field (> 80 G). [Teskel (|1974p compared power 

spectra of velocities in plages with velocity power spectra of quiet photosphere 
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and suggested that photospheric oscillations are not gravity waves. IBalthasarl 
(|1990p found that the power of five-minute range oscillations can both increase 
and decrease in faculae.The analysis of line-of-sight (LOS) velocity oscillations in 
Fei 15 648 A and 15 652 A spectral lines ( Muglach, Solanki, and Livingston, 1995 ) 
revealed dominant oscillatory signal in the velocity data is due to five-minute 
range oscillations. Also temporal variations of the magnetic-field strength were 
observed which exhibit a possible five- minute and nine - ten- minute oscillation 
( [Muglach, Solanki, and Livingston, 1995[ ). 

IKhomenko et al.\ ((MI5|) and [Centeno, Collados, and Trujillo Bueno| 
found that the time delay between photospheric (Sii 10 827 A) and chromo- 
spheric (Hei 10 830 A) five-minute oscillations in faculae is 300-500 seconds. In 
their opinion, this is a evidence of linear vertical propagation of acoustic waves. 
IKobanov et al.\ (I201ip noticed that the fivc-minutc LOS-velocity oscillations 
observed in Ha and Fe i 6569 A show ambiguous phase lags in faculae. 

This article is the third in a series dealing with the investigation of features 
of oscillations in faculae. It was preceded by [Kobanov and Pulyaev| (|2007p and 
[Kobanov and Pulyaev[ (|201ip . hereafter Paper I and Paper II, respectively. In 
this article we focus on the determination of time lag between signals of the LOS 
velocity measured in the Hei 10 830 A and Sii 10 827 A spectral lines. 

Direct measurements of the time lag between chromospheric and photospheric 
signals of the LOS velocity are of interest to test the assumption that oscillations 
propagate upwards. 



2. Instruments and Methods 

The observational spectral data analyzed here were obtained with the horizontal 
solar telescope of the Sayan Solar Observatory. The photoelectric guider of the 
telescope is capable of tracking the solar image with an accuracy of 1 arcsecond 
for several hours of observations and compensates for the rotation of the Sun 
( [Kobanov et al, 2009t . We used a Princeton Instruments RTE/CCD 256II cam- 
era 256 X 1024 pixels with a pixel size of 24 microns. One pixel corresponds to 0.24 
arcsecond along the entrance slit of the spectrograph and about 12 mA along 
the dispersion. Thus, the spectral snapshot contains information about a spatial 
region of about 60"xl". The real resolution of the telescope is 1-1.5" because 
of the Earth's atmosphere. Seeing and stray light were assessed according to the 
contrast of granulation and the limb in the visible. The intensity of stray light in 
the infrared range at the altitude of the observatory (2000 m) is low. For 10 830 A 
its intensity is roughly 16 times less than for 5000 A. We made observations in 
the Hei 10 830 A and Sii 10 827 A spectral lines. 

The LOS velocity in the Hei and Sii lines was measured by the method 
that is known as the "Doppler compensation method" ( Nikulin, Severny, and 
Stepanov, 1958). A pair of virtual slits were placed symmetrically on the spec- 
trogram on both sides of the line centre (±0.2 A for Hei and ±0.18 A for Sii) in 
such a way that the corresponding intensities Jbiuo and /red could be the same. 
Intensity in Hei was calculated as /= (/rod + /biuc)/2/c, where Ic - continuum 
intensity. When the spectral line was shifted, the virtual slits were displaced in 
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Figure 1. Increasing depth He 1 10 830 A in solar faculae (facula 5). a)Fragment of spectrogram 
and profiles of He I 10 830 A in different image parts; b)grey-scale spatial— time diagram of He I 
intensity. 

the same direction, equalizing the intensities. The displacement is proportional 
to the Doppler velocity. The initial position of the virtual slits (zero point) was 
determined relative to the telluric line H2O 10 832 A, that was used to eliminate 
the spectrograph noise. 

With a decrease in distance between the virtual slits (up to ±0.1 A), no signif- 
icant changes in He i signal were noted. Only insignificant changes in signal/noise 
ratio, associated with a decrease in the steepness of the working part of the line 
profile, were noted. For Si 1 10 827 A the picture was similar. For the both spectral 
lines behaviour of the LOS- velocity signal did not change appreciably. The width 
of virtual slits was 0.06 A for both lines. 



3. Observations 

Compared to undisturbed regions, the depth of the Hei 10 830 A line profile in 
faculae sharply increases (Figure[T^). According to our observations, this increase 
in different faculae varies from 2 to 5. This feature of the Hei line is very useful 
for precise pointing at faculae near the disk centre. At the beginning, we chose an 
object with the use of images in the CaiiH line and determined its coordinates. 
Then we slowly scanned the chosen region and corrected its position in the 
spectrograph slit, using the working camera and taking the maximum depth of 
the Hei 10 830 A line into consideration. We also should take into account the 
fact that the depth of this line increases in filaments. That is why we used Ha 
filtergrams for additional control. In most cases, the facular region was rather 
inhomogeneous and occupied the major part of the entrance aperture. Compact 
faculae were observed more seldom. In the present analysis we have used faculae 
in which the residual intensity in the core of He I 10 830 A did not exceed 0.8 for 
the space of 20". 

It is noteworthy that space-time diagrams of the Hei 10 830 A intensity can 
be used to control position of the observed object in the spectrograph slit 
throughout a time scries (Figure [T]3) . 

In the Summer 2010, we obtained 33 time series for 24 faculae near the disk 
centre; this allowed us to avoid the influence of projection effects when comparing 
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Table 1. Information on observational data, results of measurement of time lag and K. 



No. 


Date 


Time, 


Disk 


Cadence, 


Duration, 


K 


T r tj-i 
Lag ldt\, 






UT 


location 


sec 


min 




sec 


1 


17 Jul 10 


01:40 


19°N 28°W 


3 


150 


6 


78-100 


2 


04 Aug 10 


09:48 


15°N 05°W 


4 


81 


1.2 


-12-0 


3 


05 Aug 10 


00:28 


20°N 14°W 


3.5 


96 


1.5 


20-40 


4 


09 Aug 10 


00:45 


17°N 23°E 


3 


136 


1.8 


40-80 


5 


09 Aug 10 


06:10 


18°N 22°E 


4 


98 


2.9 


10-50 


6 


09 Aug 10 


09:49 


15°N 12°E 


3 


68 


1.7 


68-77 


7 


14 Aug 10 


07:29 


13°N 05°W 


3 


198 


2.2 


30-90 


8 


15 Aug 10 


04:24 


24°N 18°E 


3 


102 


3.3 


33-39 


9 


16 Aug 10 


01:58 


32°N 00°E 


3 


55 


8 


68-78 



photosphcric and chromosphcric signals. The most complete and multiple-factor 
analysis given in the article was carried out on basis of nine time series. We 
selected those series that corresponded to better seeing and minimal blurring. 
In addition, we took into account the following criteria: the remoteness from 
sunspot no less than 60", the facular size no less than 20", and homogeneity of 
facula (the residual intensity in the core of Hei 10830 A no more than 0.8), as 
well as the duration of the series no less than 50 minutes. 

Preprocessing of spectrograms consisted of standard procedures: removal of 
the "dust" effect and determination of the flat field. Using special programmes, 
we then plotted space-time diagrams of the He i 10 830 A intensity and diagrams 
showing the spatial distribution of different-frequency modes in the 1-8 mHz 
range for photospheric and chromospheric signals of the LOS velocity. With the 
use of these diagrams, we determined facular fragments in which five-minute 
oscillations are present both in the photosphere and chromosphere. We detected 
five-minute oscillations in the 2.5-4.5 mHz band. The significance level for 
the power of LOS-velocity oscillations in this band was 95%. The LOS-velocity 
signals in the facular fragments of 2" located at a distance of 5" from bound- 
aries were then subjected to frequency filtering in the 1 mHz band centred at 
3.5 mHz. Our chosen parameters of frequency filtering corresponded to those 
of |Centeno, Collados, and Trujillo Bueno| (|2009p in order to compare results 
properly. The average amplitude of the photospheric signal was then increased 
up to the average amplitude of the chromospheric signal, and the average lag 
of the chromospheric signal relative to the photospheric one was determined 
using cross-correlation. Table[T] presents obtained lag values and the amplitude 
multiplication factor of the photospheric LOS-velocity signal K= ^ch/^ph for 
each time series. The coefficient K relates only to the LOS velocities filtered in 
the 3-4 mHz band. 
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4. Analysis Results. Discussion 

Our previous articles (Paper I, Paper II, and lKobanov et ani201ip dealt mainly 
with observations in the Ha 6563 A and Pel 6569 A spectral lines. According 



to Vernazza, Avrett, and Loeser| (|198ip . the height of the Ha line formation 



is 1550-2000 km, whereas that of the Pel 6569 A line is 150-250 km ( Parnell 
and Beckers, 1969). The difference in formation heights of these lines is about 
1500 km. It is noteworthy that, according to Centeno, Collados, and Trujillo 
Bueno ipnn^ir The difference in heights in faculae for the Hei 10 830 A and Sii 
10 827 A lines is 1500 km. Spectra of chromospheric and photospheric oscillations 
of LOS velocity in faculae in IKhomenko et aQ (|2008p and Centeno, Collados, 
and Trujillo Bueno ^009p coincide down to the smallest details in the 1-8 niHz 
range. The time lag of the chromospheric signal relative to the photospheric 
one at 3.5 mHz is 300-400 seconds. Our results for the Ha and Pel 6569 A lines 
( [Kobanov et ai, 201ip contradict those of Centeno and Khomenko. Photospheric 
and chromospheric spectra differ significantly in the low-frequency (1-2 mHz) 
and high-frequency (5-7 mHz) regions. The power of chromospheric oscillations 
at 1-2 mHz frequencies obtained in Ha is very often comparable or exceeding 
the power of five- minute oscillations. Similarity between spectra is observed only 
in the central part of the range (2.5-4.5 mHz). 

But the main problem is that five-minute oscillations observed by us in the 
Ha and Fe i 6569 A lines do not show time lags of the chromospheric signal of 



300-400 seconds, unlike [Centeno, Collados, and Trujillo Bueno| p009p . In the 



majority of the phase-difference measurements, the chromospheric signal lags 
behind the photospheric signal; the lag value is however small and is rarely as 
large as 100 seconds. In some facular fragments, the chromospheric signal at 
some time intervals was slightly in advance relative to the photospheric signal. 
This should be nominally interpreted as downward propagation (Kobanov et ai, 
201lJ! 

Explanation for these contradictions is hampered by the fact that different 
pairs of spectral lines were used in the observations. 

In this article, we have the opportunity to compare results of observations, 
obtained in the same spectral lines (Hei 10830 A and Sii 10827 A). Even prelim- 
inary analysis revealed that power spectra of the LOS-velocity oscillations looked 
very similar (Pigure [2]) . This similarity implies that the frequency of oscillations 
at the height of the Hei 10 830 A line formation is identical to that at the height 
of the Sii 10 827 A line formation. According to our observations in the Hei 
10 830 A line, five-minute oscillations are well detected in the chromosphere of 
facula. Pigure [3] showing unfiltered signal of the LOS velocity confirms this. 

Amplitudes of the photospheric and chromospheric signals should be equalized 
to determine phase delay between them when they are filtered at the central fre- 
quency of 3.5 mHz. To do this we calculated the average amplitude of signals 
(for photospheric signals) and Ach (for chromospheric signals) for the entire time 
series. The amplitude of the photospheric LOS-velocity signal was multiplied by 
K = j4ch/^ph, and both signals were then superimposed in the common figure; 
this allowed us to estimate roughly the time lag between the signals. Values 
of the amplitude multiplication factor of the photospheric LOS-velocity [K] for 
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Frequency, mHz 

Figure 2. Power spectra of LOS velocity, observed in He I (thick line) and Si I (thin line), and 
normalized to their maximum values. 
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Figure 3. Unfiltered signal of chromospheric LOS velocity in He I (facula 9). 

different faculae are in the interval from 1.2 to 8 and are presented in Table [TJ 
One could suppose that this spread of K is not random and is probably connected 
with a feature of magnetic-field structure (magnetic-field strength, filling factor, 
inclination to vertical, magnetic shadow, and the presence of an admixture of 
opposite polarity). Apparently more attention should be paid to this problem in 
future. 

To determine the lag more accurately, we shifted the photosphcric signal to 
later times in step of five seconds up to a maximum of 400 seconds and found 
the cross-correlation maximum that corresponded to the average time lag [6t] in 
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Figure 4. Phase relation between photospheric (thin line) and chromospheric (thick line) 
unshiftcd signals of LOS velocity in different faculac. The signals have been filtered in a 1 mHz 
band around 3.5 niHz. Panel numbers arc in accordance with Tablc^ 
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Figure 5. Examples of LOS-velocity signals (facula 1), filtered in different frequency bands 
around 3.5 mHz. a)l mHz band; b)0.2 mHz band. 



the series. Tabled] shows values [6t] for each analysed time series. As would be 
expected, the chromospheric signal of five- minute oscillations of the LOS velocity 
lags behind the photospheric in the majority of faculae under investigation. In 
Figures U and [SJ graphs of unshifted signals of LOS velocity are given. The time 
lag differs for different series and sometimes can vary within one time series 
(Figure EJ. Ho wever, it never reaches values presented in Centeno, Collados, and 
Trujillo Bueno ipnn^ir ^ccording to our measurements, the average lag is about 
50 seconds. In TableHl the spread of lag values obtained for different parts of each 
facula is given. The spread depicts the maximum and minimum of the lag values 
calculated for several spatial elements of the facula. Only those elements were 
taken into account in which five-minute oscillations are well detected in both 
the photosphere and the chromosphere. It is worthy of note that the wavetrain 
structure of LOS velocity excludes possible error in a value divisible by period 
when making quantitative estimation of the lag {e.g. see Figure H|). We also 
tried to narrow the frequency filtering band of fivc-minute oscillations from 1 
mHz down to 0.2 mHz. No significant difference was observed afterwards, and 
the lag did not change (Figure [5]). 

We calculated the phase- difference spectra for each facula studied (Figure [5]). 
Each mark on the Figure IH means the phase difference of LOS-velocity oscilla- 
tions (He I - Si i) for each facula 2" element along the entrance slit and particular 
frequency. Each phase-difference spectrum was plotted over the whole facula 
(from 10 to 30 points). As was to be expected, phase-difference spectra seem 
ambiguous as a result of the smallness of the average lag measured. Only in four 
spectra (1,4,7,9) weak indications of propagating oscillations can be seen. We 
do not exclude the possibility that our measurements of the photospheric LOS 
velocity contains a contribution of small non-magnetic inclusions located along 
the entrance aperture. 
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Figure 6. Phase difference between photospheric and ehromospheric signals of LOS velocity 
as a function of frequency. Each mark means the phase difference of LOS- velocity oscillations 
(He I - Si l) for each facula 2" element along the entrance slit and particular frequency. Panel 
numbers are in accordance with Tabled 



Thus, our investigation of facular oscillations in the Hei 10 830 A and Sii 
10 827 A lines enables two propositions to be made. The first proposition: we 
are in full agreement with IKhomenko et a/.l(|2008[) and Centeno, CoUados. and 
Trujillo Bucno (^2009j) about the fact that spectra of the LOS- velocity oscillations 
observed in the Hei and Sii lines look very similar. The second proposition: in 
contradiction with the results obtained by these authors: the temporal lag of 
the ehromospheric signal does not exceed a period and is about 50 seconds on 
average. 

The extreme similarity between spectra can be interpreted as a result of the 
fact that all oscillations from the upper photosphere successfully reach the upper 
chromosphere. At the same time, this similarity can be caused by closer heights 
of the Hei and Sii line formation in facular regions. Possible decrease in the 
heig ht of the Hei 10 830 A line formation in faculae was noticed bv Livshits 
et al. (|1976p . In lKobanov et a/.l(|201ip we showed how similar oscillation spectra 
in faculae can be when the difference in heights is only 300-400 km, using 
observations in the Bail 4554 A and Fei 4551 A lines. Noteworthy is that the 
spatial averaging over the entire facula in this case does not make the spectra 
dissimilar. 

If the velocity of the upward wave propagation is 4-6 kms~^ (as in Centeno, 
Collados, and Trujillo Bueno J^OOBp . the difference in heights of the Hei and Sii 
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line formation in faculae will be about 250 km, considering our measured time 
lag of about 50 seconds. At the same time, if the difference in formation heights 
of these lines is 1500 km (as in |Centeno, Collados, and Trujillo Buenoll^OO^ . the 
velocity of the upward wave propagation will be about 30 kms~^. This velocity 
is six times greater than that mentioned above. 

Earlier the same problems occurred while researching oscillations observed 
at two levels of the undisturbed atmosphere. IMein and MeinI (|1976p and IMeinI 
(|1977p explained a high phase velocity of about 30 km s~^ in the chromosphere by 
restoring forces in the presence of magnetic fields. IFleck et al.\ p994p . analyzing 
observations in He 1 10 830 A and Mg i 8807 A, discussed the two alternatives that 
we cited above for faculae. In addition, IFleck and Deubnerl (jl989p considered a 
variant in which it is supposed that the measured phase-difference forms in 
the range between temperature minimum and 1000 km. They suggested that 
the o scillations become standing at heights of more than 1000 km. Gouttebroze 
et al. (|1999p found that variations of intensitv in several UV spectral lines, 
observed with the SUMER onboard SOHO, demonstrated features of evanescent 
or standing waves. 

If it is assumed that, within the cavity observed, both standing and prop- 
agating waves exist at the same time (quite an ordinary situation for any real 
resonator being the source of oscillations), then the phase lag being measured will 
depend on the contribution of each of the components. Their rate is determined 
by the degree of penetrability of the boundaries of the resonator that in real 
conditions changes following local changes in temperature, pressure and the mag- 
netic field. Previously, IKobanov et al.\ (|201ip indicated such a possibility when 
explaining the inadequacy of phase lags obtained from observations in Ha and 
Fei 6569 A. Our preliminary measurements of phase difference between velocity 
and intensity oscillations showed that the lag of 90° is more characteristic of the 
signals measured in Hci 10 830 A than for those measured in Sii 10 827 A. In 
the case of acoustic waves, this may imply that standing waves make a major 
contribution to oscillations observed at the level of He i 10 830 A formation. 

There is a contradiction that can be resolved with the supplementary informa- 
tion specifying both formation depths of these lines and a possible propagation 
velocity of oscillations. In the future it is essential to obtain useful information 
on the structure and dynamics of the magnetic field in the faculae being studied. 
It would be useful to investigate the differences in spectral oscillations of faculae 
in active regions (AR) and faculae outside AR (for example, in polar regions) 
in more detail. It is desirable also to carry out a cycle of observations of limb 
faculae simultaneously with spectroscopic and filter methods. In any case, one 
should consider the results of our analysis carried out for nine faculae. 

5. Conclusion 

We carried out an investigation of features of periodic oscillations in nine fac- 
ula regions for two heights (Hei 10 830 A and Sii 10 827 A). Power spectra of 
photospheric and chromosphcric oscillations of the LOS velocity were found to 
be very similar. According to our measurements, the lag of the chromosphcric 
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signal relative to the pliotospheric one varies from -12 to 100 seconds and is, on 
average, about 50 seconds for five-minute oscillations. 

The features that we have revealed can be interpreted in two ways: either a 
difference in heights of the He i and Si i lines formation in faculac that is substan- 
tially smaller than the generally accepted estimate (1500 km), or propagation 
velocity of fivc-minute oscillations above faculac is several times greater than 
the commonly used value of 4-6 kms~^. Another possible explanation supposes 
that within the cavity being observed the standing and propagating waves exist 
at one and the same time. In this instance the phase lag measured will depend 
on the contribution of each of the components. 

In any case, this raises new questions as to whether or not the existing models 
correspond to modern observations. 
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